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ABSTRACT 


The  results  of  a theoretical  study  of  the  interpreta- 
tion of  precursors  to  PP  (A~100°)  in  terns  of  scattering 
by  small-scale  random  irregularities  in  the  crust  and 
uppermost  mantle  are  presented.  These  results  are  compared 
with  observational  data  from  2 nuclear  explosions  recorded 
at  the  Warramunga  array  and  3 earthquakes  recorded  at  the 
NORSAR  array.  This  comparison  reveals  that  the  scattering 
interpretation  offers  a plausible  and  adequate  explanation 
of  the  detailed  observational  data,  namely ; ons?t  times, 
duration,  slowness,  azimuth  and  amplitude  variations.  The 
scattering  interpretation  of  PP  precursors,  along  with  a 
similar  interpretation  of  precursors  to  PKPPXP,  removes 
the  need  for  postulating  sharp  reflecting  discontinuities 
in  the  uppermost  few  hundred  kilometers  of  the  crust  and 
upper  mantle. 
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1 INTRODUCTION 

The  use  of  spatial  arrangements  of  seismometers  has  con- 
tributed greatly  to  detailed  investigations  of  Earth  structure 
for  more  than  a decade.  Of  particular  importance  has  been 
the  ability  to  obtain  estimates  of  slowness  (s)  and  azimuth 
( directly  from  array  data.  Such  estimates  are  inherently 
more  satisfactory  than  indirect  estimates  (e.g.,  from, 
smoothed  travel  time  observations)  and  have  been  mainly 
responsible  for  recent  advances  in  knowledge  of  the  P 
(and  S)  velocity  distributions  both  with  depth  and  laterally 
inside  the  Earth.  Direct  (s,Q)  estimates  recently  acquired 
a new  significance  with  the  advent  of  the  interpretation  of 
certain  seismic  s Lgnals  in  terms  of  scattering  by  random 
inhomogeneities.  For  scattered  signals,  slowness  (which, 
as  usual,  refers  to  slowness  or  ray  parameter  of  an  arrival 
at  the  receiver  ) will  not  generally  be  the  same  as  the 
gradient  (dT/dA)  of  an  associated  travel  time  (T, A)  curve, 
and  must  therefore  be  measured  directly.  In  addition, 
scattered  waves  will  generally  arrive  at  a receiver 
simultaneously  from  a range  of  different  directions,  and 
so  reliable  estimates  of  the  slownesses  and  azimuths  of 
the  various  contributions  must  by  necessity  involve  the 
use  of  seismic  arrays. 

In  this  paper  array  estimates  of  (s,$)  are  brought  to 
bear  in  a critical  analysis  of  an  interpretation  of  the  so- 
called  precursors  to  PP  in  terms  of  scattering  by  small-scale 
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random  inhomogenei-.ies  in  the  earth's  crust  and  upper 
mantle  recently  proposed  by  Cleary  and  Haddon  (1973)  and 
Cleary  et  al.  (1974).  Numerical  results  derived  from  the 
theoretical  treatment  of  scattering  developed  by  Haddon 
(1974)  and  Haddon  and  Cleary  (1974)  are  presented  and 
subsequently  compared  with  observational  data  from  the 
Warramunga  (WRA)  and  NORSAR  arrays.  The  NORSAR  data  and 
analysis  are  described  in  detail  because  they  have 
an  important  bearing  on  the  proper  interpretation  of  the 
entire  teleseismic  P coda,  and  no  such  detailed  results 
have  been  previously  reporced  in  the  literature.  The  com- 
parison of  theoretical  and  observational  results  presented, 
along  with  other  available  data,  confirm  that  the  precur- 
sors to  PP  (90  <A<110°)  are,  indeed,  likely  to  originate 
by  scattering  of  P waves  by  random  inhomogeneities  in  the 
crust  and  uppermost  mantle.  We  remark  that  since  the  PP 
precursors  are  a major  feature  of  the  teleseismic  P coda, 
the  results  presented  herein  provide  further  strong  support 
for  the  scattering  interpretation  of  the  entire  P coda. 


PERSPECTIVE 

The  "multiplicity  of  phases  recorded  between  P and  PP 
in  t’ e diswnce  range  95°<A<111°"  were  first  interpreted 
by  Nguyen-Hai  (19^3)  as  being  of  two  types.  The  first  type, 
whose  travel  time  curve  appeared  parallel  to  that  of  P 
itself,  was  interpreted  as  resulting  from  top-side  reflections 
from  postulated  upper  mantle  discontinuities.  This  inter- 
pretation has  subsequently  been  adopted  to  explain  various 
coda  arrivals  by  Husebye  and  Madariaga  (1970),  Davies  and 
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Frasier  (1970),  Davies  et  al.  (1971),  and  Gutowski  and 
Kanasewich  (1974)  notation  of  the  form  pdpP,  after 
Davies  et  al.  (1971)  is  now  widely  used  to  describe  the 
ray  paths  in  question.  Nguyen-Hai's  second  type  of 
coda  arrivals  were  associated  with  a travel  time  curve 
apparently  parallel  to  the  PP  curve,  and  were  interpreted 
as  resulting  from  underside  reflections  at  upper  mantle 
discontinuities.  Bolt  et  al.  (1968)  concluded  that  the 
division  of  the  observations  into  two  types  was  not  warranted^ 
they  preferred  the  interpretation  in  terms  of  reflections 
from  the  undersides  of  horizontal  discontinuities  at  depth  d, 
say, below  the  Earth's  surface,  and  adopted  the  notation  PdP. 
The  travel  time  data  presented  by  Bolt  et  al.  (1968)  and 
Bolt  (1970)  is  consistent  with  d being  generally  less  than 
400  km. 

In  subsequent  array  studies,  Wright  and  Muirhead  (1969), 
Wright  (1972)  and  Angoran  and  Davies  (1972)  have  shown, 
however,  that  the  observed  precursors  generally  have  slow- 
nesses significantly  less  than  the  valuer  entailed  by  the 
P^P  hypothesis.  Wright  and  Muirhaad  (1969)  and  Wright  (1972) 
showed  that  their  observed  travel  times  and  slownesses  are 
consistent  with  an  interpretation  involving  asymmetric 
reflection  from  postulated  dipping  crustal  interfaces  at 
distances  near  20°  from  the  source.  In  addition  to  the 
lower  slownesses,  Wright  (1972)  also  observed  slownesses 
significantly  greater  than  those  entailed  by  the  P^P  hypo- 
thesis, which  he  showed  to  be  consistent  with  a similar 
asymmetric  reflection  interpretation  involving  reflection 


4 


at  distances  near  20°  from  the  recording  station.  Although 
consistent  with  the  observed  times  and  slownesses,  the 
asymmetric  reflection  interpretation  fails  to  account  for 
the  ubiquity  of  the  precursors,  and  more  importantly  the 
duration  and  amplitude  along  observed  precursor  wave  trains. 
The  scattering  interpretation  will  be  shown  to  be  consistent 
with  the  data  in  all  of  these  respects. 


3.  THE  SCATTERING  INTERPRETATION 

The  scattering  interpretation  of  the  precursors  assumes 
that  energy  is  scattered  (i.e.,  deflected)  from  ordinary 
PP  ray  paths  by  the  presence  of  relatively  small  random 
inhomogei.eities  in  the  crust  and  upper  mantle,  with  normal 
reflection  generally  taking  place  at  the  Earth's  surface 
(see  Fig.  1) . The  amplitudes  of  the  scattered  waves  will 
depend  in  general  on  the  characteristics  of  the  random 
medium  assumed,  the  primary  wave  amplitudes  at  the  points 
of  scattering,  the  wave  frequency,  the  angular  deflections 
of  the  scattered  wave  rays,  the  volume  of  the  scattering 
region  contributing  to  the  signal  at  any  particular  time 
and  the  usual  focusing  and  defocusing  effects  associated 
with  the  Earth's  velocity  structure. 

The  scattering  and  asymmetric  reflection  interpretations 
are  similar  insofar  as  the  lower  slowness  observations  are 
associated  with  scattering  (cf.  reflection)  of  the  relatively 
high  amplitude  P waves  incident  at  the  Earth's  surface  near 
20°  from  the  source  (i.e.,  associated  with  the  so-called 
"20°  discontinuity").  The  higher  slowness  observations 
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are  associated  with  conjugate  or  dual  paths  in  which  waves 
scattered  (cf.  reflected)  at  distances  near  20°  from  the 
receiver  tend  to  be  focused  towards  the  receiver  by  the 
upper  mantle  velocity  structure.  An  essential  difference 
between  these  two  interpretations  stems  from  the  fact  that 
the  scattering  process  need  only  involve  relatively  minor 
deflections  from  normal  PP  ray  paths.  Such  deflections 
would  generally  occur  at  each  point  of  the  postulated  scatter* 
ing  region  and  a multiplic ity  of  asymmetric  ray  paths 
would  generally  contribute  towards  the  energy  arriving  at 
a receiver  at  each  particular  time,  the  totality  of  such 
contributions  giving  rise  to  the  continuous  P coda.  The 
paths  would  not  generally  lie  in  the  diametral  plane  through 
source  and  receiver  (as  represented  schematically  in  Fig.  1). 
Apparently  discrete  coda  "phases"  reported  by  previous 
authors  consist  simply  of  portions  of  the  coda  which  are 
re  .atively  <.  icrgetic  compared  with  immediately  preceding 
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portions.  On  the  scattering  interpretation  such  circumstances 
would  arise  naturally  as  a random  effect,  resulting  from 
high  constructive  interference  between  scattered  wave  con- 
tributions arriving  from  different  parts  of  the  scattering 
region,  as  well  as  from  effects  of  regional  variations  in 
upper  mantle  structure.  Strong  evidence  indicating  the 
random  nature  of  PP  precursors  is  shown  for  example  in 
figures  8 and  10  (to  be  considered  in  more  detail  later), 
where  apparently  discrete  "bursts"  of  energy  on  individual 
traces  show  little  or  no  correlation  with  similar  bursts 
on  traces  recorded  on  identical  seismometers  only  a few 
kilometers  away.  Similar  P coda  variability  across  the 
LASA  array  has  been  reported  by  Davies  (1971) . In  this 
paper  we  use  the  term  precursor  to  refer  to  the  whole  of 
the  extended  energetic  portion  of  the  P coda  preceding  PP, 
including  any  apparently  discrete  phases  which  may  occur. 


Seismic  wave  scattering  in  regions  of  small-scale 
random  heterogeneity  in  the  lowest  mantle  has  been  used 
to  explain  precursors  to  the  core  phase  PKIKP  (120°<A<143°) 
by  Haddon  (1972),  Cleary  and  Haddon  (1972),  King  et  al. 

(1973,  1974)  and  Haddon  and  Cleary  (1974).  In  close  analogy 
with  the  present  case,  conjugate  paths  exist  for  the  pre- 
cursors to  PKIKP,  and  slownesses  appropriate  to  both  paths 
have  been  observed  (Doornbos  and  Vlaar,  1973;  King  et  al,, 
1973,  1974). 

Scattering  within  the  crust  and  upper  mantle  has 
previously  been  used  to  explain  the  coda  of  local  earthquakes 
(Aki,  1969,  Aki  and  Chouet,  1974;  Dainty  and  Anderson,  1972) 
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as  well  as  short-period  P wave  time  and  amplitude  anomalies 
at  the  LASA  and  NORSAR  arrays  (Aki , 1973J  Capon,  1 74J 
Capon  and  Berteussen,  1974J  Berteussen  et  al.,  1974) t 


The  scattering  interpretation  of  PP  precursors  completely 
removes  the  need  for  postulating  sharp  reflecting  discon- 
tinuities in  the  uppermost  few  hundred  kilometers  of  the 
crust  and  upper  mantle.  In  this  respect,  the  interpretation 
differs  from  studies  too  numerous  to  mention.  Vie  remark 
that  many  of  these  previous  studies  have  implicitly  dis- 
criminated against  the  possibility  of  scattering  type 
interpretations  in  their  data  selection  and  analysis  methods. 
Indeed,  the  observational  data  to  be  presented  serves  to 
illustrate  the  difficulties  likely  to  be  encountered  in 
investigating  precursors  to  PP,  or  other  coda  "phases", 
using  data  not  sampled  densely  (and  simultaneously)  in 
space.  As  an  example,  the  recent  study  by  Gutowski  and 
Kanasewich  (1974)  involves  data  from  an  array  of  5 instru- 
ments distributed  over  an  aperture  of  almost  200  km. 

The  use  of  such  a configuration  for  studies  within  the  coda 
is  quite  unsuitable  for  the  analysis  of  short  period  waves 
having  a scattering  origin. 

4 THE  SCATTERING  MODEL  ASSUMED 

The  Earth  models  to  be  examined  are  assumed  to  be 
spherically  symmetrical  except  for  small  random  variations 
in  density  (p)  and  elastic  parameters  (k,  y)  inside  the 
uppermost  100  km  of  the  crust  and  mantle.  Models  incor- 
porating the  Jeffreys  (1939)  P and  S velocity  distributions 
ar.d  a slightly  simplified  SMAK  II  model  (Mereu  et  al.,  1974) 
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will  be  examined,  the  latter  in  order  to  demonstrate  the 
consequences  of  using  a more  detailed  representation  of 
the  upper  mantle  transition  zone.  The  random  variations 
AP,  Ak  and  Ap  in  p,  k and  u inside  the  scattering  region 
are  arbitrarily  assumed  to  have  root-mean-square  values 
of  1%  and  to  be  all  proportional  to  a common  stationary 
random  function  having  a Gaussian  autocorrelation  function 
with  correlation  distance  o=12  km.  Such  a scattering 
medium  is  broadly  consistent  with  the  conclusions  of  Aki 
(1973)  and  Capon  (1974),  who  interpreted  observed  wave- 
front  deviations  at  LASA  in  terms  of  random  variations 
in  P velocity  having  an  r.m.s.  deviation  from  the  mean 
of  about  2 to  4%  down  to  depths  of  60  to  136  km.  Comparable 
results  have  been  obtained  from  observations  of  the  kinematic 

and  dynamic  behavior  of  P waves  at  NORSAR  by  Berteussen  et  al 
(1974)  . 


Under  certain  simplifying  assumptions  (see  Haddon,  1974; 
Haddon  and  Cleary,  1974)  the  theoretical  mean-square  ampli- 
tudes of  first  order  scattered  P waves  derived  from  a simple 
harmonic  primary  P wave  by  inhomogeneities  inside  a region 
V of  an  otherwise  spherically  symmetrical  Earth  model  are 
given  approximately  by 


<e2>  = 


K 
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where  k - 2tt/ A denotes  the  wavenumber  of  the  primary 
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and  all  the  other  symbols  are  as  defined  in  Haddon  and 
Cleary  (1974).  In  deriving  the  numerical  results, 
the  primary  wave  envelope  function  E(aQT)  has 
been  assumed  to  have  the  simple  form 

0 t <0 

1 0<x<10<sec 

0 t> 10  sec 

corresponding  to  a simple  harmonic  wavelet  cf  duration  10 
seconds.  Further  details  of  the  scattering  theory  are 
given  in  Haddon  and  Cleary  (1974) . 


E (aQx  ) = 


5 THEORETICAL  RESULTS 

On  the  random  medium  hypothesis,  scattered  waves  will 
generally  arrive  simultaneously  at  a receiver  R from  a wide 
range  of  different  directions.  Fig.  2 shows,  for  example, 
the  approximate  theoretical  ranges  of  slowness  and  azimuthal 
direction  for  scattered  waves  of  PP  type  arriving  at  a 
receiver  at  an  epicentral  distance  A=100°  at  the  times  indi- 
cated for  the  Jeffreys  velocity  model.  The  scattered  wavelets 
, arriving  at  R at  any  particular  time  are  seen  to  consist  of 
two  groups,  one  having  slownesses  in  the  range  4.5  to  7.5 
sec/deg  corresponding  to  scattering  at  points  nearer  to 
the  source  than  to  the  receiver,  and  the  other  having  slow- 
nesses in  the  range  7.5  to  13.5  sec/deg  to  scattering  at 
points  nearer  to  the  receiver  than  to  the  source  (see  Fig.  1) . 
The  theoretical  slownesses  of  these  two  groups  will  be  com- 
pared with  observational  data  in  a subsequent  section. 

The  two  sets  of  curves  in  Fig.  2 are  conjugate  to  one 
another  in  the  sense  that  scattered  wave  rays  corresponding 
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to  points  on  any  curve  in  the  lower  part  of  the  figure  are 
associated  with  primary  wave  rays  corresponding  to  points 
on  the  curve  for  the  same  time  in  the  upper  part  of  the 
figure  and  vice  versa.  It  may  he  noted  that  the  in 

the  curves  between  slownesses  of  4.5  and  5.1  sec/deg  in  the 
lower  part  of  the  figure  are  associated  with  the  cusps  in 
the  Jef freys-Bullen  P travel  time  curves  near  20°. 


Fig.  2 Ranges  of  slowness  (s/deg)  and  azimuthal  deviation 
(deg)  of  waves  scattered  near  the  surface  of  the  earth  and 
subsequently  arriving  at  a receiver  at  A=100°  at  the  times  indicated 
for  a Jeffreys  velocity  model.  The  upper  ( >7 . 5 s/deg)  and 
lower  (■'7.5  s/deg)  siowness  groups  correspond  to  scattering 
on  the  receiver  and  source  sides  of  the  paths  respectively. 

Note  the  wide  azimuthal  spread  of  wave  scattered  on  the 
receiver  side  of  the  symmetrical  reflection  point. 
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Fig.  3 shows  the  corresponding  relative  theoretical 
r.m.s.  amplitudes  of  the  scattered  waves  arriving  at 
A=95  , 100°  and  105°  for  waves  of  period  1.0,  1.5,  2.0  and 
3.0  sec  as  derived  by  numerically  integrating  Eq.  (1)  for 
the  assumed  model.  The  relative  amplitude  of  the  associated 
direct  primary  PP  arrival  as  derived  using  simple  ray  theory 
(Bullen  Ch.  8)  is  approximately  2.0  units.  It  must  be  em- 
phasized vhat  although  the  theoretical  r.m.s.  amplitude 
curves  show  smooth  amplitude  variation,  random  fluctuations 
in  amplitude  along  individual  wave  trains  are  of  course  en- 
tailed by  the  scattering  meet  m assumed. 

The  relatively  sharp  increase  in  scattered  wave  ampli- 
tudes for  A=105°  near  time  t=1035  sec  shown  in  Fig.  3 cor- 
responds to  the  arrival  of  waves  scattered  from  the  relatively 
high  amplitude  P waves  incident  at  the  Earth's  surface  at 
distances  near  20°  from  the  source  (i.e.,  from  near  the 
so-called  20°  discontinuity)  and  from  waves  following  con- 
jugate paths  where  scattering  takes  place  about  20°  from 
the  receiver.  Following  the  amplitude  increases  from  t=1035 
sec  to  about  1050  sec,  the  theoretical  r.m.s.  scattered  wave 
amplitudes  remain  moderately  large  right  up  until  the  ar- 
rival time  of  the  main  PP  phase.  Similar  comments  apply 
to  the  results  for  the  other  two  distances  given.  The 
resuits  given  have  been  obtained  assuming  a surface  focus,* 
the  results  would  not,  however,  be  expected  to  be  significantly 
different  for  focal  depths  of  up  to  100  km  or  so.  For 
much  deeper  foci  the  results  may  differ  significantly, 
because  for  foci  at  depths  >413  km  in  the  Jeffreys 
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Fig.  3 Theoretical  r.m.s.  amplitude  variations  along  scatt-'" 

inS  at/=2i * * * 5  ' 100  and  105  £or  the  Jeffreys  verity 

modei  and  simple  harmonic  primary  wave  trains  of  period  l.o/ 

an^  l:°  se?s  and  10  sec  duration.  The  associated 

ampiitude  of  the  primary  PP  phase  as  derived  using  simple 

ray  theory,  is  2.0  units.  y H 
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model,  for  example,  the  primary  wave  amplitude  at  the  surface 
of  the  model  would  be  quite  different  from  that  for  a sur- 
face focus,  since  the  T-A  curve  is  single-valued  in  the 
former  case. 

The  effects  on  the  theoretical  results  of  replacing 
the  Jeffreys  upper  mantle  velocity  distribution  by  a sub- 
stantially different  upper  mantle  distribution  as  repre- 
sented, for  example,  by  a s lightly  simplified  version  of  the 
SMAK  II  model  of  Mereu  et  al.  (1974)  are  now  examined,  since 
it  is  now  known  that  the  velocity  gradients  in  the  upper 
mantle  transition  zone  can  differ  significantly  from  those 
in  the  Jeffreys  distribution.  The  Jeffreys  and  the  simpli- 
fied SMAK  II  velocity  distribution  used  in  the  calculations 
arr  plotted  for  comparison  in  Fig.  4.  The  amplitude  and 
(s,$)  calculations  involved  in  Figures  2 and  3 were  re- 
peated for  the  SMAK  model,  and  the  corresponding  results 
are  presented  as  Figures  5 and  6. 

An  important  difference  between  the  results  for  the 
two  models  arises  from  the  shift  of  some  4°  to  shorter 
distance  of  the  caustic  point  C associated  with  the  "400  km 
discontinuity"  in  the  SMAK  model  (see  Fig.  4)  due  mainly 
to  the  shallower  depth  (315  km)  of  this  discontinuity  as 
compared  to  the  Jeffreys  model.  This  results  in  the  build- 
up of  amplitude  in  the  theoretical  precursor  wave  trains 
occurring  up  to  about  20  seconds  earlier  for  the  SMAK  model 
than  for  the  Jeffreys  model  (cf.  Figs.  3 and  6).  Another 
important  difference  (cf.  Figs.  2 and  5)  is  that  in  place 
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Fig.  4 Comparison  of  the  Jeffreys  (heavy  line)  and  a 
slightly  simplified  SMAK  II  (lighter  line)  upper  mantle  P 
velocity  models  used  in  theoretical  calculations.  Note  that 
the  caustic  cusps  C in  the  associated  reduced  travel  time 
d^f£^s  occur  near  A=19  for  the  Jeffreys  model  but  near 
A=15  for  the  SMAK  II  model. 


of  the  single  group  of  higher  slowness  scattered  wave  arrivals 
for  the  Jeffreys  model  there  are  two  distinct  groups  in  the 
SMAK  II  model,  corresponding  to  the  two  sets  of  T-A  branches 
associated  with  the  abrupt  changes  in  velocity  gradient  near 
315  and  670  km  in  the  latter  (Fig.  4) . The  differences  be- 
tween the  (s,$)  results  shown  serve  to  indicate  the  extent 
to  which  regional  variations  in  upper  mantle  structure  may 
affect  the  slownesses  and  azimuths  of  scattered  wave  arrivals. 
In  particular,  the  differences  indicate  that  the  uncertain- 
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ties  do  not  warrant  a very  precise  interpretation  of  slow- 
ness observations  in  the  higher  slowness  range  (>7.5  sec/ 
deg),  as  will  become  apparent  later.  The  theoretical 
results  given  here  will  be  compared  with  observational 
data  in  a subsequent  section. 


T«IO(mc) 


1000  (M 

TRAVEL  TIME  (secs) 


Fig.  6 Theoretical 
wave  trains  for  A=95 
model  (cf.  Fig.  3). 


Qr .m. s^  amplitude  variations  along  scattered 
, 100  and  105  for  the  SMAK  II  velocity 
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DATA  AND  ANALYSIS 

Th„  analysis  techniques  used  in  the  compilation  of 
vatijnai  data  suit-able  for  a neaningful  comparison  with 
theoretical  results  are  crucially  important  in  the  present 
context,  and  so  will  be  described  in  detail  below.  sin<  e 
the  analysis  of  the  WRA  data  to  be  examined  has  been  dis- 
cussed eisewhere  (Cleary  et  al.,  1974),  the  present  section 

will  oe  confined  to  the  selection  and  analysis  of  NORSAR 
data. 


6.1  Data  S • lection 

The  phenomenon  under  examination  manifests  itself  over 
■nany  seconds  and  wide  (<,.♦)  ranges  (e.g.,  see  Figures  2 and 
5).  It  is  essential  that  for  each  event  examined,  the 
processing  "space"  is  of  comparable  dimensions  and  is 
sampled  in  adequate  detail.  These  requirements  present 
us  with  a computational  problem  of  unusual  magnitude,  and 
so  we  have  restricted  attention  to  3 events.  Since  the 
full-aperture  NORSAR  array  begain  to  operate  on  a routine 
basis  in  1971,  events  of  suitable  magnitude  and  in  a suit- 
able distance  range  for  analysis  (mm>5.6,  90°<A<110°)  have 
occurred  in  the  two  distinct  source  regions  framed  in 
Fig.  7.  One  event  in  the  Northern  Chile  region  (in  region 
A,  Fig.  7)  and  two  events  from  the  Molucca  Passage  (in 
region  B,  Fig.  7)  were  selected  for  analysis,  the  hypo- 
central  information  for  which  is  listed  in  Table  1.  These 
events  allow  us  to  examine  in  detail  precursor  wave  trains 
from  two  widely  separated  source  regions,  and  also  to 
make  comparisons  between  data  from  within  a relatively 
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before  PP  than  have  some  channels  with  much  lower  P ampli- 
tudes (02B,  03B) . The  small  power  level  variations,  as 
well  as  the  extended  duration  of  the  PP  precursor  wave  train, 
suggests  immediately  that  the  precursor  generating  mechanism  is 
effective  over  a large  region  of  the  Earth's  upper  mantle, (i .e. , 
a large  area  of  slowness-azimuth  space).  All  these  features,  as 
well  a-  the  arrival  time  of  the  energy  relative  to  the  P arrival 
time,  cannot  be  satisfactorily  explained  by  structural 
complexities  in  the  immediate  vicinity  of  NORSAR,  but  are 
fully  consistent  with  the  expectations  of  the  scattering 
interpretation.  Also  consistent  with  theoretical  expecta- 
tions (briefly  referred  to  in  an  earlier  section)  is  the 
relative  lack  of  coherence  between  even  closely  separated 
subarrays  (see  NORSAR  configuration,  inset  in  Fig.  7) . In 
fact  the  recordings  from  the  6 instruments  within  individual 
7-kn  subarrays,  such  as  those  plotted  in  Fig.  9,  display 
only  quasi-coherence,  and  are  suggestive  of  extensive  inter- 
ference between  a plethora  of  arrivals  arriving  from  dif- 
ferent directions.  Note  that  the  traces  in  Fig.  9 are  each 

normalized  independently,  and  have  been  filtered  in  the 
frequency  range  0.6  Hz  to  3.0  Hz. 

Fig.  1 0 is  comprised  of  selected  data  from  subarray 
central  instruments  for  events  1 and  2 (Table  1) . In  the 
time  interval  of  particular  interest  (i.e.,  two  or  so 
minutes  preceding  PP) , although  the  traces  in  this  figure 
differ  in  detail  both  within  the  array  aperture  for  a single 
event  and  also  between  the  separate  events,  all  of  the  data 
presented  appears  to  exhibit  the  same  general  character. 
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Fig  9.  Caption  at  bottom  of  following  page 


We  wxll  demonstrate  shortly,  however,  than  in  respect  of 
the  (s,4>)  structure,  the  data  from  the  two  different  regions 
are  very  different.  Fig.  10  includes  one  trace  (labelled 
LG)  plotted  more  as  it  might  appear  on  the  recording  from 
a conventional  "analogue"  station.  This  particular  trace 
is  dominated  by  the  strong  core  phase  arrival  (PKiKP) , and 
the  details  of  the  coda  and  PP  itself  are  lost.  This  example 
serves  to  illustrate  the  importance  in  studies  of  this  kind 
of  having  magnetically  recorded  data  and  the  associated 
capability  to  increase  arbitrarily  the  gain  in  plotting. 

6 . 2 Data  Analysis 

The  processing  procedure  adopted  for  investigation 
of  the  distribution  of  power  in  time,  slowness  and  azimuth 
involved  vhe  formation  of  some  588  array  beams  for  each  second 
of  data  analysed.  The  search  space  adopted  was  4 to  12  sec/deg 
for  slowness  and  $expected  + 60°  for  azimuth,  although  the 
slowness  range  12  to  15  sec/deg  was  also  investigated  on  an  ex- 
ploratory basis  for  selected  time  intervals.  Intra-subarray  time 
delays  were  calculated  and  subarray  beams  formed  for  12 
(s,4>)  reference"  directions  distributed  symmetrically 
within  this  search  space  by  summing  the  6 appropriately 
phased  traces  for  each  of  the  22  subarrays.  Using  the 
subarray  beams  for  each  reference  direction  in  tur"  the 
inter-subarray  delays  were  varied  as  appropriate  for  direc- 


Fig.  9 Precursor  section  of  seismograms  from  ail  6 of  the 
short  period  instruments  of  the  7-km  aperture  subarrays  01A 
and  01B  for  event  3.  The  plot  start  time  is  1 hr  0 mins. 

Each  of  the  traces  shown  has  been  filtered  in  the  range  0.6  Hz 
to  3.0  Hz  and  normalized  independently  of  the  others. 
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tions  corresponding  to  each  point  in  turn  of  a symmetrical 
grid  of  7 x 7 (s,<t>)  values  centred  on  each  reference  direc- 
tion. The  beam  "power"  (energy  per  second)  was  evaluated 
for  each  grid  point  by  summing  the  phased  subarray  beams 
to  form  a full  array  beam,  and  then  summing  squared  full 
array  beam  amplitudes  over  20  samples  (1  second) . This 
procedure  yielded  beam  power  at  each  point  of  a 14  x 42  grid, 
with  spacing  0.62  sec/deg  and  3.1°  in  s and  <t>  respectively 
for  each  second  of  data  analysed.  With  regard  to  systematic 
errors  in  the  (s,<t>)  values,  the  time  delays  appropriate  to  the 
(s,4)  values  defininq  the  grid  points  were  "decalibrateu" 
in  the  processing  using  the  NORSAR  empirical  time  delay 
correction  base  (Berteussen,  1974) . Since  this  correction 
base  is  near  optimal,  the  decalibration  procedure  should 
effectively  remove  any  serious  distortion  in  the  grid  point 
locations  arising  from  structure  underlying  NORSAR. 


The  processing  procedure  described  was  used  to  process 
the  105  consecutive  seconds  of  events  2 and  3 immediately 
preceding  and  including  the  PP  arrival,  and  the  120  seconds 
of  event  1 preceding  the  PKiKP  arrival.  In  addition,  a 
five  second  noise  sample  preceding  P,  P itself  and  some 
85  seconds  of  the  "early"  ccda  of  event  3 were  also  analysed 
for  comparison  with  the  precursor  section  of  coda.  Since 
a single  15  second  data  interval  required  a minimum  of  90 
minutes  processing  time  on  the  IBM  360/40  computers  at  the 
NORSAR  Data  Processing  Center,  the  enormity  of  the  computa- 
tion problem  requires  no  further  emphasis.  Alternative 
processing  procedures  were,  of  course,  considered,  but  lacked 
sufficient  robustness  in  the  complex  signal  situation  at 
hand. 
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The  values  of  power  at  each  (s,4>)  grid  point  for  each 
second  of  the  wave  train  were  outputed  on  cards  or  magnetic 
tape  for  subsequent  contouring  in  a separate  processing 
step.  In  order  to  facilitate  the  interpretation  of  power 
peaks  for  subsequent  comparison  with  theoretical  results 
such  as  those  presented  in  Fig.  2,  contour  plotting  in  the 
slowness-azimuth  plane  of  power  value  averaged  over  2 or 
3 seconds  was  found  to  be  most  suitable  for  digestion  of 
all  results.  A typical  slowness-azimuth  power  plot  for  a 
3 second  data  segment  of  total  travel  time  1025  secs  (Event  2) 
is  illustrated  in  Fig.  11.  Vespa-type  plots  (Davies  et  al., 
1971)  in  the  diametral  azimuth,  but  with  some  control 
exercised  ov?r  the  array  azimuthal  response  by  incorporating 
data  from  adjacent  azimuth-time  planes  in  an  averaging  pro- 
cedure, were  also  plotted  (see  Fig.  16) . 

6 • 3 Interpretation  of  Slowness-Azimuth  Diagrams 

The  processing  described  effectively  transforms  the 
array  data  into  a convenient  representation  in  slowness- 
azimuth  space  by  the  application  of  a family  of  matched 
spatial  filters.  It  is  essential  to  give  consideration  to 
the  response  of  these  filters  in  slowness,  azimuth  and  time 
in  order  that  any  information  extracted  from  the  transformed 
data  can  be  reliably  ascribed  as  a "true”  manifestation  of 
Earth  structure  rather  than  as  an  artefact  of  the  processing. 
We  are  concerned  in  particular  with  establishing  criteria 
for  assigning  significance  to  the  various  peaks  in  the  con- 
tour plots  of  power  in  the  slowness-azimuth  planes  (e.g.. 

Fig.  11)  . 
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Fig.  11  Contoured  beam  power  levels  in  slowness-azimuth 
space  averaged  over  3 seconds  for  a mean  travel  time  of 

J2Hns®c. f0rVKVenJ  2;  The  contours  represent  power  levels 
in  dB  below  the  absolute  maximum.  The  3 dB  contour  is  darkened 

of  P°wer  levels  between  2 and  1 dB  and  1 and  0 dB 
are  hatched  and  dotted  respectively.  Theoretical  slowness- 

curves  for  scattered  wave  arrivals  for  a travel  time 
for  A=10°  for  both  the  Jeffreys  (solid  line) 
and  SMAK  II  (broken  line)  are  plotted  for  comparison  with 
the  observed  power  di  -ibution.  Note  that  the  locations 

P°wef  1 shown  are  in  close  agreement 
with  the  theoretical  curves.  Also  note  the  wide  spread 
of  azimuths  of  observed  power  peaks  in  this  example.  This 
hKS  ^een/frec3uently  observed  in  the  present  study 

(cf.  Figs  il"  l^anHsK*  l0Wer>  Sl°WneSS  9rOUP  o£  ^rivals 
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The  beamforming  response  of  the  NORSAR  array  to  a mono- 
chromatic steady— state  signal  has  been  discussed  by  several 
authors  and  is  conveniently  represented  in  wavenumber  space 
(see  inset  to  Fig.  7).  The  response  to  broadband  transient 
signals  is  a smoothed  version  of  the  steady-state  response 
(Kelly,  1967,  Steinberg,  1971),  but  can  differ  also  in  the 
mean  sidtlobe  level  (see  Steinberg,  1971).  In  a typical 
case  observed  at  NORSAR,  the  v/orst  sidelobe  in  the  response 
of  a bean  power  analysis  on  the  main  P signal  from  Event  2 
(Table  1)  was  measured  to  be  slightly  less  than  10  dB  below 
the  main  lobe  level,*  this  is  some  5 dB  worse  than  the  level 
indicated  by  the  steady-state  theory.  Such  discrepancies 
are  probably  due  mainly  to  the  variations  in  signal  ampli- 
tude and  shape  within  the  array  aperture,  which  violate 
certain  fundamental  assumptions  in  the  steady-state  calcula- 
tions. Indeed,  Doornbos  and  Husebye  (1972)  have  demonstrated 
that  amplitude  variations  between  individual  sensors  can 
have  the  effect  of  introducing  spurious  side  lobe  peaks 
into  the  array  response.  Such  effects  are  expected  to  be 
troublesome  in  the  analysis  of  PP  precursor  wave  trains 
insofar  as  the  individual  contributions  to  the  wave  train 
are  subjected  to  similar  distortions  as  the  main  P signal 
by  structure  beneath  the  array. 

The  response  of  the  array  to  a composite  signal,  in- 
volving several  interfering  contributions  of  various  ampli- 
tudes and  arrival  directions , cannot  be  simply  estimated 
from  a superposition  of  appropriate  steady-state  responses 
and  is  extremely  complex.  In  addition,  the  time  response 
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of  the  array,  which  depends  on  the  duration  and  pulse  shape 
of  a transient  arrival,  will  act  to  restrict  the  resolution 
of  interfering  components.  Doornbos  and  Husebye  (1972) 
investigated  some  of  these  effects  using  synthetic  data 
and  elucidated  the  difficulties  in  interpreting  arn.y  beam 
information  in  the  presence  of  multiple  arrivals.  In 
slowness-azimuth  space  the  situation  is  further  complicated 
by  the  fact  that  the  response  pattern  for  each  arrival 
separately  depends  on  the  particular  arrival  direction. 

In  view  of  these  difficulties,  we  examined  our  data 
directly  in  order  to  establish  empirically  an  "interpreta- 
tion threshold"  in  the  contoured  power  plots  of  2 to  3 
second  data  intervals.  This  examination  revealed  that  inter 
pretation  of  power  peaks  more  than  2 to  3 dB  below  the  most 
energetic  power  peak  would  be  problematic  in  that  it  would 
greatly  increase  the  number  of  peaks  per  time  interval. 

In  view  of  the  array  time  response  limitations,  it  seems 
unreasonable  to  extract  multiple  peaks  within  a 2 to  3 
second  interval,  except  perhaps  in  the  situation  where 
peaks  of  very  similar  power  are  present,  and  so  this 
threshold  is  a purely  practical  one.  We  remark  that  if 
we  chose  to  interpret  all  peaks  down  to  6 dB,  say,  many 
of  the  measured  values  would  then  find  no  explanation  in 
the  scattering  or  any  other  suggested  interpretation. 
However,  at  this  level  all  the  "spurious"  peaks,  without 
exception,  occur  in  time  intervals  containing  a more 
powerful  "valid"  peak,  which  suggests  that  the  numerous 
peaks  3 to  6 dB  below  the  main  peak  are  likely  to  be 
artefacts  of  the  array  response  in  an  extremely  complex 
signal  situation. 
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It  is  relevant  to  comment.  on  the  distinction  between 
our  interpretation  threshold  as  compared  to  that  suggested 
by  Capon  and  Davies  (1971),  who  investigated  the  specific 
problem  of  interpreting  contoured  beam  power  diagrams 

produced  with  a processor  similar  to  that  used  in  the 

. 

present  study.  These  authors  utilized  techniques  of 
statistical  pattern  recognition  in  establishing  criteria 
for  detecting  significant  arrivals  (i.e.,  power  peaks) 
within  a signal-generated  coda,  and  concluded  that  a 
discrete  seismic  arrival"  could  be  declared  present  for 
a minimum  contrast  of  9 dB  between  the  main  peak  power 
and  the  average  background  noise  level  just  prior  to  signal 
arrival.  Although  these  authors  used  this  criterion  to 
examine  arrivals  within  the  teleseismic  P coda,  we  differ 
insofar  as  their  assumed  threshold  :.s  based  on  the  dis- 
tinction between  so-called  signal  and  supposedly  independent 
preceding  noise.  Such  a distinction  is  inappropriate  in  the 
present  investigation  in  which  we  are  seeking  to  determine 
approach  directions  of  contributions  having  relatively  high 
intensity  within  wave  trains  which  consist  of  sums  of  con- 
tributions of  similar  intensities  arriving  from  widely  dif- 
ferent directions  (cf . Figs.  2 and  5)  . In  such  circumstances 
all  the  contributions  arriving  from  directions  within  the 
main  lobe  of  the  array  response  pattern  appropriate  to  a 
Par^icular  direction  may  be  thought  of  as  constituting  the 
"signal"  corresponding  to  that  direction.  All  contributions 
not  encompassed  by  the  main  response  lobe  then  effectively 
constitute  the  "noise".  Since  the  main  response  lobe  cor- 
responds to  a range  of  directions  much  smaller  than  the 
total  range  of  directions  from  which  the  various  contribu- 
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tions  could  arrive  (cf.  Figs.  2 and  5),  the *teignal-to-noise 
ratio"  is  likely  to  be  generally  poor  and  so  well-defined 
high  intensity  peaks  would  not  be  expected.  It  i_s  to  be 
expected,  however,  that  absolute  maximum  energy  peaks  ob- 
tained in  a search  of  the  entire  (s,$)  space  would  usually 
correspond  to  directions  of  relatively  high  scattered  wave 
intensity,  although  exceptions  may  sometimes  occur  because 
of  some  of  the  effects  discussed  above. 

In  summary,  beam  power  values  have  been  evaluated 
at  each  point  of  an  extensive  grid  of  (s,$)  values  for 
each  of  many  consecutive  seconds  of  data  from  three  events. 
These  data  have  subsequently  been  combined  over  2 and  some- 
times 3 seconds,  and  power  values  contoured  in  the  slowness- 
azimuth  plane.  The  largest  power  values  in  continguous 
15  second  time  intervals  were  used  as  reference  levels  in 
the  appropriate  consecutive  plots.  The  slowness  and  azi- 
muth values  corresponding  to  no  more  than  two  significant 
peaks  in  each  contour  plot  were  recorded  and  will  be  com- 
pared with  the  theoretical  expectations  of  the  scattering 
interpretations  in  the  following  section. 


7.  COMPARISON  OF  OBSERVATIONS  WITH  THEORY 

In  this  section,  observational  data  from  the  WRA  and 
NORSAR  arrays  will  be  compared  with  appropriate  theoretical 
results.  We  have  chosen  to  use  theoretical  amplitude 
results  for  o=12  km,  T=3.0  secs  (cf.  Figs.  3 and  6)  for 
comparison  with  observed  amplitude  variations  despite  the 
fact  that  the  dominant  periods  in  the  observational  data 
are  significantly  less  than  3 secs  (<■'0.3-1. 7 secs).  We 
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do  not  consider  this  apparent  mis-match  to  be  serious, 
however,  because  the  relative  amplitudes  of  the  observed 
PP  precursors  are  much  too  large  for  the  linear  scattering 
theory  to  be  strictly  applicable.  Indeed,  the  dominance  of 
precursor  (scattered)  energy  over  direct  PP  energy  in  our 
observational  data  (see  Figs.  8,  10,  12)  indicates  that 
scattering  effects  are  so  strong  that  there  can  be  no  dis- 
tinction between  primary  and  scattered  waves.  In  these 
circumstances  results  from  the  linear  theory  may  still  be 
used  to  model  consequences  of  the  postulated  scattering 
process,  but  the  parameters  a and  T do  not  then  have  the'* 
usual  physical  significance,*  the  applicability  of  linear 
scattering  theory  in  such  circumstances  is  discussed  more 
fully  in  a later  section. 

I 

7 . 1 Observations  at  Warramunga 

Fig.  12  shows  processed  traces  of  two  Novaya  Zemlya 
explosions  recorded  at  WRA  (A  = 106°)  on  Oc  c’-sr  27th,  1966, 
and  October  14th,  1970,  compared  with  corresponding 
theoretical  r.m.s.  amplitude  curves  for  scattered  waves 
derived  for  a Jeffreys  model  and  a period  of  3.0  seconds. 
Note  that  i)  the  amplitudes  of  the  observed  PP  precursor 
wave  trains  start  to  increase  to  relatively  significant 
values  about  85  secs  before  the  arrival  of  the  main  PP 
phase,  and  ii)  the  observed  amplitudes  remain  significant 
right  up  until  the  arrival  of  PP.  In  both  of  these  respects 
the  observational  data  are  in  good  agreement  with  the 
theoretical  amplitude  variation  curves. 


Fig.  12  Processsed  WRA  array  traces  showing  PKiKP , PP  and 
precursors  to  PP  from  Novaya  Zemlya  explosions  (A=106°)  on 
27  October  1966  (example  (a))  and  14  October  1970  (example 
(b) ) . The  traces  shown  are  (from  bottom)  time  channel, 
single  seismometer  output,  delayed  summation  of  both  red  and 
blue  line  seismometers  (ZR,ZB) , array  beam  (ZR+ZB)  and  one- 
second  time  averaged  product  (SIG  ✓ (ZRxZB) ) . Sections  (a) 
and  (b)  have  been  aligned  according  to  travel  time.  Note 
the  duration  of  the  precursor  wave  trains,  and  also  the  marked 
difference  in  the  arrival  time  of  precursor  energy  maxima 
from  tgo  explosions  which  had  almost  identical  location 
(A=106  ).  Theoretical  r.m.s.  amplitudes  for  A=105  , T=3.0 
sec,  for  the  Jeffreys  velocity  model  are  plotted  for  com- 
parison with  the  data. 


Comparison  of  the  two  observed  wave  trains  shown  in 
Fig.  12  also  provides  clear  evidence  of  the  random  nature 
of  precursor  wave  trains,  as  was  originally  pointed  out 
by  Cleary  et  al.  (1974) . The  notable  dissimilarities  in 
the  wave  trains  from  two  sources  in  a similar  location  are 
extremely  difficult  to  account  for  on  any  previously  sug- 
gested interpretation  of  the  precursors,  but  are  fully 
consistent  with  the  random  variations  to  be  expected  with 
the  scattering  mechanism  when  the  source  locations  are 
separated  by  a distance  comparable  with  or  larger  than  the 
characteristic  size  of  the  inhomogeneities  present. 
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bleary  et  al.  (1974)  have  shown  that  the  slowness  o£ 
the  precursor  wave  train  for  the  1970  explosion  migrates 
to  higher  values  with  Increasing  time.  The  corrected 
slowness  values  (see  Cleary  et  al.,  1974)  for  the  signals 
some  80  and  40  seconds  before  the  PP  arrival  are  plotted  on 
the  appropriate  theoretical  slowness-azimuth  diagram  in  Fig 
13.  Also  plotted  on  this  figure  are  some  results  for  an 
event  at  a distance  of  105.5°  from  the  Yeiio„knlfe  array 


presented  by  Wright  (1972),  these  results  representing  the 
first  reported  precursor  slownesses  significantly  higher  than 
the  PP  slowness.  The  agreement  between  observed  and  theoretical 
results  in  Fig.  13  is  fully  satisfactory  in  all  cases. 
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M?^nred  WRA  (s'*>  values  (crosses)  for  the 

i thCth^"  H197°TWre^ain  (example  (fa)'  gig-  12)  compared 
m ShithSu  fc  1 (s,<f)  results  f°r  A=1056  for  the  Jeffreys 
model.  The  measured  slowness  values  are  labelled  with  their 
corresponding  measured  travel  times.  Also  shown  (dots  in 
lrcies)  are  some  slowness  measurements  made  at  the  Yellow- 
knife array  by  Wright  (1972)  for  his  event  11. 


37 


7 • 2 Observations  at  NORSAR 
Event  1 


Six  subarray  central  instrument  traces  and  their  cor- 
responding short-term- average  (STA)  traces,  obtained  by 
simple  rectification  and  smoothing,  are  plotted  in  comparison 
with  the  theoretical  r.m.s  amplitude  curves  for  both  the 
Jeffreys  and  SMAK  models  in  Fig.  14.  The  theoretical  curves 
correspond  well  with  the  general  features  of  the  observational 
data.  No  comparison  with  the  PP  arrival  is  possible  as  a 
result  of  the  strong  PKiKP  arrival  (see  Fig.  10).  The 
observed  travel  tines  in  Fig.  14  have  been  corrected  for 
focal  depth  (~105  km)  by  the  addition  of  15  seconds,  which 
is  half  the  measured  P-pP  separation  (see  Fig.  10a) . The 
slowness-azimuth  measurements  corresponding  to  the  corrected 
travel  times  are  plotted  in  Fig.  15.  The  high  slowness 
range  of  the  theoretical  results  has  been  excluded  sinc° 
no  significant  energy  peaks  having  slownesses  in  this  range 
were  found.  The  grouping  of  measurements  into  20  second 
time  windows  reflects  the  fact  that  a more  precise  interpre- 
tation is  not  warranted. 

Scrutiny  of  Fig.  15  shows  that  the  (s,$)  data  is  in 
satisfactory  agreement  with  the  theoretical  results.  In 
particular,  the  small  observed  deviations  from  "true" 
azimuth  correlate  well  with  the  lower  slowness  values  ob- 
served, this  effect  is  remarkably  consistent  with  theoretical 
expectations  as  indicated  by  the  concave  downward  curvature 
of  the  (s,$)  curves  for  the  lower  slowness  group  in  Fig.  2. 
Since  all  measured  slownesses  are  within  + 15°  of  the 
diametral  azimuth,  the  observational  data  can  be  conveniently 
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Fig.  14  Seismogram  traces  (filtered  in  the  range  0.6  to 
3.0  Hz)  and  associated  short  term  average  (STA)  traces, 
showing  PP  precursor  wave  train  at  6 selected  NORSAR  sub- 
array central  instruments  (as  labelled)  for  event  1 in 
comparison  with  relative  theoretical  r.m.s.  amplitudes  for 
the  Jeffrgys  (broken  line)  and  SMAK  II  (full  line)  models 
^or  and  T=3.0  secs  (see  Figs.  3 and  6)  . The  amplitude 

scale  for  the  theoretical  results  has  been  chosen  arbitrarily 
for  satisfactory  correspondence  with  observed  amplitudes,  but 
is  the  same  in  all  parts  of  the  figure. 
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respondin^  to  each  of  these  time  intervals  are  hatched  on 
the  appropriate  part  of  the  figure.  natened  on 
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summarized  in  a Vespa-type  plot.  Such  a plot  is  presented 
as  Fig.  16,  but  this  plot  differs  from  a conventional 
Vespagram  insofar  as  the  data  from  6 adjacent  azimuth  ranges 
has  been  averaged  in  the  calculation  of  the  power  values. 

The  migration  of  values  to  higher  slowness  with  increasing 
tine  is  clear.  The  bounds  of  the  lower  slowness  group  in 
the  theoretical  results  included  in  Fig.  15  are  superimposed 
on  Fig.  16  to  facilitate  comparison.  The  power  values  at 
slownesses  >8  sec/deg  are  almost  exclusively  8 dB  or  more 
below  the  peak  powers  at  lower  slownesses . 

Events  in  Begion  B 

In  Fig.  17  the  theoretical  r.m.s.  amplitudes  for  both 
the  Jeffreys  and  SMAK  models  are  compared  with  the  frequency 
filtered  and  STA  traces  from  6 subarray  central  instruments 
for  each  of  Events  2 (Fig.  17a)  and  3 (Fig.  17b) . The  ob- 
servational data  agrees  satisfactorily  with  both  theoretical 
curves,  although  the  gradual  energy  build  up  after  about 
970  seconds  is  perhaps  in  better  agreement  with  the  results 
derived  for  the  SMAK  model. 

The  measured  (s,4>)  values  for  both  of  these  events 
are  plotted  on  theoretical  curves  for  the  SflAK  model  in 
Fig.  18a  to  e.  The  SMAK  model  is  preferred  for  this 
comparison  since  it  proved  marginally  superior  in  explain- 
ing the  observational  data.  Both  models  accomodate  the 
observational  data  reasonably  well,  however,  and  a discrimi- 
nation between  various  upper  mantle  velocity  distributions 
on  the  basis  of  our  data  is  not  really  justified.  In  all, 

76  data  points  are  plotted  in  Fig.  13;  more  than  one  half 
of  these  are  in  complete  agreement  with  theoretical  expecta- 
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Fig.  16  Vespa-type  plot  of  the  disti 
function  of  slowness  and  travel  time 
gram  has  been  formed  by  averaging  beai 
adjacent  azimuths  (in  the  ra 
and  each 
azimuth  of  246  , 
depth,  r J 
hatched  (2  dB 
to  the  power  maxima 
the  figure.  The  bounds  of  s 
group  of  arrivals  for  A=105 
are  plotted  (heavy  lines) 

values.  Note  that  the  i 

theoretical  upper  bound  line 


secong  of  data.  The  data  plotted  is  ior  a nominal 
anc*  the  time  axis  is  corrected  for  focal 
Major  powet  Peaks  are  dotted  (1  dB  <pover  < dB)  and 
" -3  "Power  <1  dB),  the  power  values  being  relative 

in  the  intervals  indicated  at  the  top  of 
fewness  for  the  lower  slowness 
for  the  Jeffreys  velocity  model 
for  comparison  with  the  observed 
measured  slownesses  closely  follow  the 
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Fig.  18  Measured  (s,4>)  values  at  NORSAR  for  event  2 (triangles) 
and  event  3 (crosses)  plotted  in  comparison  with  theoretical 
(s,<D)  values  for  A=100  for  the  SMAK  II  model.  The  measured 
values  are  grouped  in  depth-corrected  travel  time  intervals 
of  20  secs,  and  the  theoretical  (s,<J>)  results  appropriate 
to  each  of  the  time  intervals  are  hatched  on  the  appropriate 
frames  of  the  figure.  Note  the  wide  spread  of  measured  azimuths 
associated  with  higher  slowness  (s>7.5  s/deg)  observations 
(cf . Fig.  11) . 


tions,  and  only  7 of  the  remaining  observed  values  are  in 
serious  conflict  with  the  theoretical  results.  In  view  of 
the  considerable  difficulties  in  making  (s,4>)  measurements 
within  the  P coda,  as  evidenced  by  the  paucity  of 
published  data,  v:e  regard  the  observational  data  as  being 
in  very  satisfactory  agreement  with  the  theoretical  expecta- 
tions. Random  interference  effects  between  the  various 
contributions  together  with  distortion  of  the  signal  by 
structure  beneath  the  array  may  account  for  the  inconsistent 
observations,  although  we  have  not  investigated  this  pos- 
sibility  in  detail.  Additional  features  of  the  data  presented 
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which  support  the  scattering  interpretation  of  the  PP  pre- 
cursor wave  trains  are  the  increase  in  the  spread  of  the 
azim  ithal  observations  associated  with  the  higher  slowness 
observations,  and  the  gradual  coalescing  of  the  high  and 
low  slowness  regimes  with  increasing  tine.  The  extent  of 
the  spread  of  the  power  in  azimuth  is  also  apparent  in  Fig.  11, 
as  is  the  consistency  with  the  theoretical  curves  included  in 
this  figure.  Finally,  we  remark  that  the  very  low  coherence 
in  the  highly  energetic  PP  precursor  wave  trains  analysed 
here  is  in  itself  quite  significant  support  for  the  scatter- 
ing interpretation. 

DISCUSSION 

Theoretical  calculations  for  the  assumed  models  indicate 
that  on  the  linear  scattering  theory  the  total  energies  in  the 
"higher"  and  "lower"  slowness  groups  of  scattered  wave  arrivals 
at  a receiver  are  practically  equal.  However  the  ranges  of  slow- 
ness and  azimuth  for  the  lower  slowness  group  are  very  much 
smaller  than  for  the  higher  slowness  group  (see  Figs.  2 and  5), 
and  so  the  lower  slowness  group  would  be  much  more 
coherent  at  a receiver  array.  Thus,  other  things  being 
equal,  the  highest  scattered  wave  intensities  would  be 
expected  to  be  associated  with  the  lower  slowness  arrivals. 
These  theoretical  expectations  are  borne  out  well  by  the 
WRA  data  and  the  NORSAR  data  from  region  A (Fig.  7) . In 
addition,  while  lower  slowness  PP  precursors  have  been 
commonly  observed  (see,  for  example,  references  in  introduc- 
tion), only  Wright  (1972)  has  previously  reported  observa- 
tions of  the  higher  slowness  group.  We  suggest  that  ob- 
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servations  of  the  higher  slowness  group,  such  as  the  exten- 
sive NORSAR  observations  from  region  B (Fig.  7),  are  most 
likely  to  be  made  when  the  lower  slowness  group  is  rela- 
tively weak  or  absent.  This  could  happen  as  a result  of 
large-scale  regional  variations  in  the  upper  mantle  scatter- 
ing structure,  and  comparison  of  NORSAR  (s,<i>)  data  from 
regions  A and  B is  strong  evidence  for  such  regional  varia- 
tions. There  is  also  a suggestion  in  the  (s,<J>)  data  from 
region  B (see  Fig.  18)  that  scattering  on  the  western 
side  of  the  diametral  azimuth  is  somewhat  stronger  than 
on  the  eastern  side.  We  interpret  the  dominance  of  the 
higher  slowness  group  in  the  NORSAR  data  from  region  B as 
offering  strong  evidence  for  the  presence  of  widespread 
random  irregularities  in  the  vicinity  of  the  Ural  mountains, 
some  20  to  30  from  NORSAR.  This  conclusion  is  broadly 
consistent  with  implications  of  the  WRA  data,  which  pro- 
vides independent  evidence  for  effective  scattering  in 
a similar  region. 

In  the  observed  wave  trains  shown  in  Figures  12,  14 
and  17,  the  mean  amplitudes  of  the  PP  precursors  are  approx- 
imately equal  to  the  amplitudes  of  the  main  PP  phases.  The 
corresponding  theoretical  precursor  amplitudes  for  the 
model  are  only  about  a quarter  of  the  amplitude  of  the 
PP  phase.  Increasing  the  assumed  proportional  r.m.s.  varia- 
tions in  density  and  elastic  parameters  from  1%  would  in- 
crease the  theoretical  r.m.s.  amplitudes  of  the  scattered 
waves  relative  to  the  main  PP  (primary)  phase  proportionally. 
It  is  evident  from  the  theoretical  amplitude  curves  shown, 
however,  that  even  for  1%  variations  a significant  part  of 
the  energy  in  the  primary  P wave  train  would  be  transformed 
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scattered  waves  in  the  assumed  model  (particularly  for 
shorter  period  waves) , an  effect  not  allowed  for  on 
linear  scattering  theory,  while  this  indicates  a certain 
inadequacy  of  the  linear  theory  in  the  present  context,  scrutiny 
of  the  theory  shows  that  Eq.  (1)  may  still  be  applied  to 
model  the  scattering  process  if  it  is  assumed  that  the  energy 
in  each  narrow  beam  of  primary  wave  rays  entering  the  (thin) 
scattering  layer  is  scattered  so  that  the  resultant  relative 
directional  intensity  of  the  energy  emerging  from  the  layer 
may  be  approximated  empirically  by  exp  f-a2sin2 4>/2 ] , where 
* denotes  scattering  angle  and  a is  a constant.  Further  assuming 
ordinary  ray  theory  to  give  the  relative  amplitudes  of  the 
primary  waves  entering  the  postulated  scattering  layer  and 
the  effects  of  geometrical  spreading  on  the  scattered  waves 
after  emerging  from  the  layer  (as  in  the  linear  theory),  Eq.  (1) 
would  give  the  relative  r.m.s.  amplitudes  of  the  resultant 


scattered  waves  of  PP  type  arriving  at  a surface  receiver  in 
the  assumed  model.  In  these  circumstances  the  kinematic 
characteristics  of  the  scattered  waves  (i.e.,  slownesses, 
azimuths  and  arrival  times)  would  be  the  same  as  for  the 
linear  theory,  but  the  parameters  T , k and  a appearing  in 
Eq.  (1)  would  not  have  any  physical  significance.  The  use 
of  Eq.  (1)  to  represent  strong  scattering  effects  thus  amounts 
essentially  to  postulating  a plausible  directional  intensity 
function  for  scattered  waves  emerging  from  each  element  or 
the  postulated  scattering  layer  rather  than  postulating  the 
physical  characteristics  of  the  layer  and  deriving  the  con- 
sequent scattering  effects.  While  in  the  former  case  there 
is  no  established  connection  between  the  characteristics 
of  the  postulated  random  medium  and  the  consequent  scattering 
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effects  produced,  this  in  no  way  affects  our  basic  interpre- 
tation of  the  PP  precursors.  Notwithstanding  all  the  various 
uncertainties  involved,  the  extent  of  agreement  between  the 
observed  and  theoretical  results  presented  strongly  supports 
the  scattering  interpretation  of  PP  precursors. 

9.  CONCLUSIONS 

Seismic  scattering  by  random  inhomogeneities  in  the 
crust  and  upper  mantle  has  been  shown  to  account  adequately 
for  detailed  observational  data  on  precursors  to  PP  from 
the  WRA  and  NORSAR  arrays.  These  data  included  principally 
onset  times  and  duration,  slrwness,  azimuth  and  amplitude 
variations  of  precursor  wave  trains.  Many  of  these  features 
j-ind  no  explanation  in  any  previously  suggested  interpreta- 
tion of  the  precursors.  A similar  scattering  interpretation 
accounts  for  precursors  to  PKPPKP  having  lead  Mnes  less 
than  about  90  seconds  (King  and  Cleary,  1974;  Haddon  et  al., 
1975) . Scattering  d also  appear  to  account  adequately 
for  the  entire  hort  period  teleseismic  P coda  (Cleary  et  al., 
1974) . With  regard  to  the  latter,  the  study  of  precursors 
to  PP  is  important  insofar  as  there  seems  little  prospect 
of  establishing  unambiguously  the  relative  contributions 
of  various  coda-generating  mechanisms  (e.g.,  Key,  1967,  1968; 
Husebye  and  Madariaga,  1970;  Davies  et  al.,  1971 J Greenfield, 
1971,  Phinney  et  al. , 1974)  from  studies  of  the  coda  between 
P and  any  precursor  arrivals.  The  scattering  interpretation 
of  PP  (and  PKPPKP)  precursors  entirely  removes  the  need  for 
postulating  sharp  reflecting  discontinuities  in  the  upper- 
most few  hundred  kilometers  of  the  crust  and  upper  mantle. 
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